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Abstract 

We investigate the linear classical stability of Bogomornyi-Prasad-Sommerfield (BPS) on three 
domain wall solutions in a system of three coupled real scalar fields, for a general positive potential 
with a square form. From a field theoretic superpotential evaluated on the domain states, the con- 
nection between the supersymmetric quantum mechanics involving three-component eigenfunctions 
and the stability equation associated with three classical configurations is elaborated. 
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I. INTRODUCTION 

The kink solution of a field theory is an example of a soliton in 1+1 dimensions 
. From field theoretic superpotential [6( evaluated on the domain states, the algebraic 
ramework of supersymmetry in quantum mechanics (SUSY QM), as formulated by Witten 
3], may be elaborated. The SUSY QM generalization of the harmonic oscillator raising and 
lowering operators has several applications [8]. The generalization of SUSY QM for the case 
of matrix superpotential, is well known in the literature for a long time. See, for example, 
for one-dimensional systems the works in Ref . [9 1 about non-relativistic quantum systems, 
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id . 1 1 11 ] about one single field and in [12l . 1 131 Il4l ] with two-field and in [15| on a three-field 



potential models in 1+1 dimensions. 

The classical configurations with domain wall solutions are bidimensional structures in 
3+1 dimensions. They are static, non-singular, classically stable Bogomol'nyi [16| and 



Prasad-Sommerfield 



17l | (BPS) soliton (defect) configurations, with finite localized energy 



associated with a real scalar field potential model. Domain walls have several applications 
in condensed matter and cosmology. The BPS states are classical configurations that satisfy 
the first order differential equations and the second order differential equations (equations 
of motion). 

Recently, marginal stability and the metamorphosis of BPS states have been investigated 
. via SUSY QM, with a detailed analysis for a 2-dimensional N = 2— Wess-Zumino 
model in terms of two chiral superfields, and composite dyons in iV = 2-supersymmetric 
gauge theories. Domain walls have been recently exploited in a context that stresses their 
connection with BPS-bound states. While Rajaraman has applied the trial orbit method 
for the equation of motion, here we use the trial orbit method for the first order differential 
equations associated to three real scalar fields. However, for solitons of three coupled scalar 
fields there are no general rules for finding analytic solutions since the nonlinearity in the 
potential increases the difficulties to obtain the solutions of the BPS equations and field 
equations. 

This paper is organized as follows: In Section II, we investigate domain walls configura- 
tions for three coupled scalar fields, and supersymmetric non-relativistic quantum mechanics 
with three-component wave functions is implemented. In Section III, a scalar potential model 
of three coupled scalar fields is investigated. Our Conclusions are presented in Section IV. 
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II. THREE COUPLED SCALAR FIELDS 



We consider the classical soliton solutions of three coupled real scalar fields in 1+1 di- 
mensions. They are static, nonsingular, classically stable and finite localized energy solutions 
of the field equations. 

Recently, a superfield formulation of the central charge anomaly in quantum corrections 
to soliton solutions with N=l SUSY has been investigated [19]. The superaction in terms 
of the superfield and the superpotential in superspace Q Q ), a = 1, 2, is written as 



S N=1 = X - J dx 2 d 2 & (j2 D&iD&i + W($!, $ 2 , $ 3 )) , 



(1) 



where D is the supercovariant derivative given by 



D = d e + i&fd„ (2) 

and the 7^ are represented in terms of the Pauli matrices, 7 = a y and 7 1 = ia x . The 
superpotential W($i, $2, $3), yields the component-field potential V(0j), where $j are chiral 
superfields which, in terms of bosonic (tpi), fermionic (ipi) and auxiliary fields (Fi), are 
G— expanded as shown below: 

$i = h + + -^Fi, (3) 

where Gj and = G* are Grassmannian variables. 

The Lagrangian density for such nonlinear system in the natural system of units (c = 
h = 1), in (l+l)dimensional space-time with Lorentz invariance, in terms of the bosonic 
fields only, is written as 

£ Ms) = \il (<V>;) 2 - Vfa), (4) 

Z 3=1 

where = = (t,z), <pj = (j)j(t,z), (j = 1,2,3) are real scalar fields and 

_ (_|- ) — ) i s the metric tensor. The potential V(4>j) = V((f>i, 02,03) is any positive 
definite functional of <f>j. The general classical configurations obey the equation bellow: 

Q2 Q2 Q 



dt 2Yj dz 2Y3 
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which, for static soliton solutions, is equivalent to the following system of nonlinear second 
order differential equations 







V, (j = 1,2,3), 



(6) 



where prime represents differentiation with respect to space variable. There is, in literature, a 
trial orbit method for finding static solutions for certain positive potentials, which constitutes 
a "trial and error" technique 

We can get the masses of the bosonic particles, using the results above, from the second 
derivatives of the potential: 



ml = 



d 2 V 



\z— >±oo 



% = 1,,2,3. 



(7) 



A. Linear Stability and SUSY QM 



Since the potential V(0j) is positive it can be written in the following square form, 
analogous to the case with only a single field 6], 



Vfa) = Vfa, 02, 03) = - E ^(01> 02, 03), Ujtfx, 02, 03) 

1 j=l 



dW 



, (j = 1,2,3). (8) 



In this case, one can write the total energy 



E 



+oo I 

dz- 

oc 2 



^l) 2 + (02) 2 + (03) 2 + 2K(0 4 ) 



in the BPS form, consisting of a sum of squares and surface terms, 



E 



00 / 1 1 1 8 N 

} dz -(0; - Uxf + -(0; - U 2 f + i(0' 3 - f/ 3 ) 2 + |V 



(9) 



(10) 



which are always positive. Thus, the lower bound of the energy (or classical mass) is given 
by the fourth term, viz., 



E > 



'°°dz^-w[M^,Mz),Uz)] 

X) OZ 



(11) 



where the superpotential W = W[4>i(z), 02 (2), 03 (2)] wu ^ be discussed below. The BPS 
mass bound of the energy resulting in a topological charge is given by 
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E BPS = T ij = \W[M j ]-W[M i ]\, 



(12) 



where M, and Mj represent the BPS vacuum states and are the extrema of W.. 

In this case the BPS states satisfy the following set of first order differential equations 
associated to three real scalar fields: 



$3 = 02,03) • 



(13) 



Now, let us analyze the classical stability of the soliton solutions in this nonlinear system, 
which can be investigated by considering small perturbations around 4>i(z), <$>i{z) and 4>z{z): 

4>j(t, z) = 4>j(z) + Vj (t, z), (j = 1, 2, 3), (14) 
where the fluctuations r/j(t, z) can be expanded in terms of the normal modes as 



Vj(t, Z) = J2 ^j,rJ]j,n{z)e W ^ n \ U ltTt = ^2,n = ^3,n = W n , 



(15) 



with ej tn being real constant coefficients. Thus, the stability equation for the field becomes 
a Schrodinger-like equation of a three-component eigenfunction \l/ n , 



where 



m n = uj 2 J> n , n = 0,l,2,---, 



n = -i^ + v F (z), 



(16) 



(17) 



with I being the 3x3 identity matrix. The term Vf(z) is the 3x3 fluctuation Hessian matrix 
and the excited mode is given by 



/ 



Vl,n( Z ) 



\ 



V V3,n( z ) J 



(18) 



Since 



v 92 v - 32 v - V 

Ftj d(pid(pj d(f)jd(pi Fjt 



(19) 
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we see that 7i is Hermitian. Hence the eigenvalues uj 2 of 7i are real. 

The Schrodinger-like equation ffl6|) and the Hessian matrix Vf(z) are obtained by con- 
sidering a Taylor expansion of the potential V((j)j) in terms of 7)j(t, z) and retaining the first 
order terms in the equations of motion. 

From (1171) we find a bilinear form of TC given by 



where 



n = a + a~ 



(20) 



A ± = ±I-^+W(z), A + = (A~)\ W\z)=W(z) 
dz 



(21) 



Using these first order differential operators that appear in analysis of classical stability 



associated to a single field 



6], we find 



A + A' 



d 2 d 2 
+ —V, 



(22) 



'33 dz 2 

which are exactly the diagonal element of 7i. Therefore, it is easy to show that the linear 
stability is satisfied, which means that 



--< U >=< A + A- >= (A-^ n ) j (A-^ n ) > 0, 



as was stated. Also, the ground state of 7i becomes 



(23) 



( UMz)]^ 

U 2 [<P(z)] (24) 

which represents the bosonic three-component zero mode. 

The 3x3-matrix superpotential satisfies the following Ricatti equation associated to the 
non-diagonal fluctuation Hessian Vp(z): 



( 



W 2 + W = V F (z) 



V Fn (z) V F12 (z) V F13 (z) y 

V F12 (Z) V F2 2{z) V F2 z{z) 

\ V F13 (z) V F2 3{z) V Fm {z) 



(25) 
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The Ricatti equation (125]) holds the BPS states only. According to the Witten' SUSY 



model 
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we have 



V + 1 W 1X6 

where ^± (z) are three-component eigenfuctions. In this case, the graded Lie algebra of the 
supersymmetry in quantum mechanics for the BPS states may be readily realized as 

■ A + A- \ (U- . 

H SUSY = [Q-,Q + ] + =\ = > ( 27 ) 

o ^-^+ y 6xg w H 

Q±L = o = (<5-) 2 = (Q+) 2 , (28) 
where Q± are the 6x6 supercharges of Witten SUSY model, viz. 

/ 



Q_ = o- ® A , Q+ = Ql 



CT-j 



(29) 



A^ 
V 

with the intertwining operators, ^4 ± , in terms of 3x3-matrix superpotential, are given by Eq. 
( ]2"6"]) and <x± = |(eri ± 202), where <ti and o% are the Pauli matrices. Of course, the bosonic 
sector of H S usy is exactly the fluctuating operator given by H~ = H = — I-j^ + \p(z), 
where V_ = Vp(z) is the non-diagonal fluctuation Hessian. The supersymmetric fluctuation 
operator partner of 7i- is 

H + = A-A + = A + A- + [A~, A + ] =H-- W'(^), (30) 

so that the SUSY partner is given by V + = V_ — ~W'(z). 
Starting with 

H ^n) = E (n)y(n) ^ A + A~^^ = E H ^ ] (31) 

and multiplying fl3T|) from the left by A~ we obtain 

A-A + (A-^ y ) = E^(A-^ y ) H + (A-^ ] ) = E^\a-^ ] ). (32) 
Since A~^^ = 0, comparison of (152"]) with 

H + ¥? = A-A+¥? = E^¥^\ (33) 



7 



leads to the immediate mapping: 



E[ n) = y(n) ^ ^("+1)^ = ^ ^ ^ ^ 

Repeating the procedure but starting with (I3"3"j) and multiplying the same from the left 
by A + leads to 

•A+zTC/^J ) = E^iA+vP), (35) 
so that it follows from (13~I]) . (151) and (j33|l that 

^L n+1) ocA+^ n) ,n = 0,1,2,.... (36) 

The intertwining operator A~(A + ) converts an eigenfunction of 7i_(7i + ) into an eigen- 
function of TL + (TLJ) with the same energy and simultaneously destroys (creates) a node of 
. These operations just express the content of the SUSY operations effected by 



Q + and Q_ connecting the bosonic and fermionic sectors of the SUSY fluctuation operator 

(EZD- 

The SUSY analysis presented above in fact enables the generation of a hierarchy of 
Hamiltonians with the eigenvalues and the eigenfunctions of the different members of the 
hierarchy in a simple manner. Calling 7i_ as 7i\ and 7i + as Ti.2, and suitably changing the 
subscript qualifications, and by repetition of the above procedure leads to the generation of 
a hierarchy of Hamiltonians given by 

U n = + V n (x) = A+A- + E$» = A-_.AU + E { X (37) 

A± = ±I-^ + W n (z), Wt{z)=W n {z) (38) 
whose spectra satisfy the conditions 

E?- 1 = E^ 2 = ... = El?\ n = 2, 3, . . . , M, (39) 
* n 1 - 1 «AtAt...AU* i n ) - (40) 



S 



Note that the nth-member of the hierarchy has the same eigenvalue spectrum as the first 
member 7ii except for the missing of the first (n— 1) eigenvalues of TCi. The (n-l)th-excited 
state of Hi is degenerate with the ground state of TL n and can be constructed with the use 
of (JHU) that involves the knowledge of Ai(i — 1, 2, . . . , n — 1) and \E^ ) . 

III. A POTENTIAL MODEL OF THREE SCALAR FIELDS 

Let us consider the following generalized potential in terms of bosonic fields only 



where a > and $ > 0. Note that the symmetry Z 2 xZ 2 is only preserved if <j) 2 = or 
Pi — P2 — condition is satisfied. When 3 = this potential becomes the two-field 
potential model recently investigated 

From ([6]) and (HT1) the equations of motion under static limit, for each bosonic component- 
field, become 




+ -eg (afa - 2/3 2 2 - a^f 



(41) 




(42) 




(43) 
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+ 2^ 2 03 [Ml - «0i0 2 - ft) 

+ </) 3 (-2/? 2 2 + OUfo. ~ OiPlf . 



(44) 



The corresponding superpotential model in field theory is given by 

W(<j>i) = ^01 - ^0? - - f 0103 + &0203 - &02 + ^103, (45) 

where a > and $ > 0. Thus, it is required that all <pj, j = 1, 2, 3 satisfy the BPS state 
conditions: 



h =-1^01+202+203- — 



> 2 = Ml - "0201 - P2 



' =-0 3 («0i-2/? 2 02 -«A). (46) 



Note that the BPS states saturate the lower bound so that Ebps = |Wy| is the central 
charge of the realization of N — 1 SUSY in 1+1 dimensions. Thus, the vacua are determined 
by the extrema of the superpotential 

dW 

—-=0, j = 1,2,3 (47) 

which provides the possible vacuum states. 

Let us now consider a projection on the (4>i,<t>2) plane in order to find an explicit form 
of domain walls. In this case, if we choose 3 = and (3 2 = 0, the superpotential W((f>i), 
becomes 

W(<fH) = ^0i - ^0? - f 0i0 2 2 , (48) 

which has been discussed recently and some orbits for this projection have been investigated 
in references Q, Q. Indeed, using the vacuum states, the following trial orbit, 



G(cf) 1 ,4> 2 ) = c 1 4>l + c 2 (f> 2 2 + c 3 = 0, (49) 

from |2 = Jg^ + Jg^/ - = &b an d the BPS states (HE}, we obtain the constant 
coefficients, q. Therefore, the ground state for this projection on the (0i, 4> 2 ) plane becomes 
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-afa{z)<k{z) (50) 

V I 

which represents the bosonic three-component zero mode. However, we have seen that 
if is a normalizable three-component eigenstate, one cannot write in terms of 
\]/^ in a similar manner to ordinary supersymmetric quantum mechanics. Also, A~ = 
ip^ -^o), is only valid for the unidimensional case with one-component eigenstate. 

There the superpotential is given by 

W(x) = -^-\nijj {0 \x). (51) 
ax 

This can not be seen in the example treated here of the classical stability analysis for three 
coupled real scalar fields. 

The matrix superpotentials in SUSY QM for other projections of the system considered 
here will be investigated in a forthcoming paper, which can not be written as Eq. ( JBTT) . 

IV. CONCLUSION 

In this paper, we consider the classical stability analysis for BPS domain walls associ- 
ated with a potential model of three coupled real scalar fields, which holds for non-ordinary 
supersymmetry (SUSY). The approach of effective quantum mechanics provides a realization 
of the SUSY algebra in the three-domain wall sector of the non-relativistic formalism. 

The tensions can be deduced from the charge central properties in a model that present 
N = 1 SUSY, which depend on the manifold of vacuum states, = |TU[Mj] — W[Mj]|, 
where M, and Mj represent the vacuum states. 

We have shown that the positive potentials with a square form lead to three-component 
non-negative normal modes u n 2 > 0, analogous to the case with a single field [6|, so that 
the linear stability of the Schrodinger-like equations is ensured. 

A general three-component zero-mode eigenfunction is deduced. The Ricatti equation 
( 1251) holds the BPS states only. A detailed analysis in SUSY QM for such of a potential 
model will be reported separately. 

In a three-field potential model there is a static classical configuration inside a topological 
soliton. We point out that the superpotential model investigated here can be applied in order 
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to implement new string junctions 20] by extendeding BPS domain walls and string theory 



of three bosonic moduli 
in 



211 ] . The set of potential BPS junctions that have been identified 



221 ] contain the junctions that appear in [23] . Also, the BPS saturated objects with 
axial symmetry (wall junctions, vortices), in generalized Wess-Zumino models, have been 
investigated in Ref. [241 ]. 
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